Quantum limit for avian magnetoreception: How sensitive can a chemical compass be? 
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The chemical compass model, based on radical pair reactions, is a fascinating idea to explain 
avian magnetoreception. At present, questions concerning the key ingredients responsible for the 
high sensitivity of a chemical compass and the possible role of quantum coherence and decoherence 
remain unsolved. Here, we investigate the optimized hyperfine coupling for a chemical compass in 
order to achieve the best magnetic field sensitivity. We show that its magnetic sensitivity limit 
can be further extended by simple quantum control and may benefit from additional decoherence. 
With this, we clearly demonstrate how quantum coherence can be exploited in the functioning of a 
chemical compass. The present results also provide new routes towards the design of a biomimetic 
weak magnetic field sensor. 



Introduction. — Quantum effects in biology has re- 
cently raised considerable interest in the communities of 
physics, biology and chemistry [IHZ]- Current research is 
concerned with understanding how long-lived quantum 
coherence may exist and what role it may play in sev- 
eral important biological processes [5HT5], The develop- 
ment of quantum biology may in turn help us to advance 
biomimetic quantum technologies. One particular exam- 
ple of quantum biology, in which quantum effects may 
be relevant, is avian magnetoreception. As one of the 
main hypotheses to explain the magnetic sensing of many 
species |16l 117] , the radical pair mechanism involves un- 
paired electrons that have intermediate correlated spin 
states [2H6]. Several important experiments have pro- 
vided compelling evidence to support the hypothesis that 
the radical pair mechanism is exploited in avian magne- 
toreception [TH] [TH] (although other mechanisms are be- 
lieved to make contributions as well), and important ob- 
servations towards the understanding of the physiological 
basis for avian magnetoreception J2HH2I]. (a) It was ar- 
gued that quantum entanglement (i.e. non-classical cor- 
relations between two radicals) [52] and interaction with 
environment may be important for the functioning of a 
chemical compass |13H15j . Inspired by the experimen- 
tal results, some design principles have been proposed 
for a chemical compass, e.g. freeing one radical pair of 
its nuclear spin environment (called reference-and-probe 
model) [231 [21], applying gradient fields [25] or by ex- 
ploiting quantum criticality [26 . It is however unknown 
what is the sensitivity limit that a chemical compass can 
achieve, whether and how quantum coherence can be ex- 
ploited to achieve maximal magnetic sensitivity. 

In this paper, we first investigate what kind of hyper- 
fine coupling in a chemical compass can lead to better 
magnetic field sensitivity. Motivated by our numerical 
optimization results, we focus on the simplest model con- 
figuration with one nuclear spin, which appears to offer 
the best sensitivity. This simple model is an approxima- 
tion of multiple nuclei environments with one effective 
dominant coupling component. It does also allow us to 
obtain analytic expressions for the optimized magnetic 



sensitivity so that we can better understand the underly- 
ing physics. A strong anisotropic hyperfine coupling has 
usually been considered to give a high sensitivity. Inter- 
estingly, we find a very rich structure in the sensitivity as 
a function of the hyperfine coupling, which has not been 
expected before. We show that a chemical compass with 
a reasonable radical pair lifetime can in principle work for 
a magnetic field as weak as a few tens of nano Tesla (nT) . 
Furthermore, we demonstrate that the sensitivity limit 
can be extended with simple quantum control, i.e. using 
an extra time-dependent magnetic field to control the 
radical pair dynamical processes. This result illustrates 
how quantum coherence can be exploited to improve the 
sensitivity of a chemical compass. Decoherence, coming 
from the interaction with the environment (out of the 
core chemical system of the radical pair and the adja- 
cent nuclei), is unavoidable at room temperature and is 
generally regarded to exert a destructive influence on a 
chemical compass. Perhaps surprisingly, we demonstrate 
that the performance of a chemical compass can be very 
robust against dephasing and can even be improved by 
increasing the dephasing rate. This phenomenon clearly 
shows that an interplay between coherent spin dynamics 
of radical pairs and environmental decoherence can play 
an interesting role in a chemical compass. 

Optimization of a chemical compass. — In a chemical 
compass, each radical has an unpaired electron coupled 
to an external magnetic field B and a few nuclei via the 
Hamiltonian [2S] H = -^ e B{S A + S D ) + Ek=A,D Ej & ' 
Tkj ■ Ik •) where 7 e = — .g e /is is the electron gyromagnetic 
ratio, TL denote the hyperfine coupling tensors, and Sk 
and Ify are the electron and nuclear spin operators, re- 
spectively. For simplicity, we assume that dipole-dipole 
and exchange interactions between two radicals can be 
neglected. This approximation is valid when the radical- 
radical distance is relatively large or the dipole-dipole 
and exchange interactions partially cancel each other, as 
proposed in Ref. |27j . The hyperfine coupling is the most 
essential ingredient of a chemical compass, as it leads to 
the coherent interconversion between the electron spin 
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singlet and triplet states, which then recombine into dif- 
ferent chemical products at a certain rate. It has not 
been fully understood what properties it should have to 
enable a chemical compass to have maximal magnetic 
sensitivity. 

We consider a radical pair created in a spin-correlated 
electronic singlet state \S) = - lit)) with the 

same singlet and triplet recombination rates, i.e. kg = 
&t = k. Under ambient conditions, the nuclear spins 
initially are assumed to be at thermal equilibrium, i.e. 
described by a density matrix as ■ Ij / dj , where dj is 
the dimension of the jth nuclear spin and lj is the iden- 
tity matrix. The singlet yield calculated from the con- 
ventional phcnomenological approach |28j is quite similar 
to the one from the quantum-measurement-based master 
equation [T5], but requires much less computational ef- 
forts. Explicitly, it is calculated as <&s = Jo°° r (t)fs{t)dt, 
where r(t) = kexp(-kt) is the radical recombination 
probability distribution, and fs(t) is the singlet ratio 
of the electron spin state p s (t) at time t, i.e. fs(t) — 
(S\ Ps (t)\S). 

As a known design principle of a chemical compass, 
partially inspired by the behavior experiments with Eu- 
rope robins [TS], it was proposed that in an optimally 
designed chemical compass, only one of the radicals has 
strong and anisotropic hyperfine interactions, while the 
other is free from the coupling with nuclear spins [231 121] • 
Thus, we take this kind of reference-and-probe model as 
the starting point for our optimization procedure. The 
target function is the magnetic sensitivity of the singlet 
yield quantified by 
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where $^ max ) and $^ min -' are the maximum and min- 
imum singlet yield for different directions of the mag- 
netic field. In this scenario our numerical optimization 
of the hyperfine couplings suggests that a chemical com- 
pass having more nuclear spins will not perform better 
than having only one nuclear spin |30) . We thus consider 
the simple case of one nuclear spin with the anisotropic 
hyperfine coupling tensor as T = diag{0, 0, a} that can 
give the best sensitivity. Such a simple model allows 
us to obtain analytic results which are quite useful for 
understanding the essential mechanism responsible for a 
chemical compass. The real molecule for avian magne- 
toreception would be more complicated than this model, 
while the physics revealed by the simple model may still 
provide insights into avian magnetoreception. In fact, it 
can be an approximation of more general models with 
one dominant hyperfine component. It may also be pos- 
sible that the motion/disorder of protein may average the 
complicated interactions, and may eventually result in a 
simple effective Hamiltonian [29] . 

The direction of the weak magnetic field B is char- 
acterized by (6 \<j>). Without loss of generality, we can 
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FIG. 1. (Color online) The magnetic sensitivity of the sin- 
glet yield, Ds, as a function of the ratio between hyperfine 
coupling a and magnetic field strength B. The radical pair 
recombination rate is k = 0.5/iS -1 , and the magnetic field is 
set as geomagnetic field B = 46/iT. A chemical compass can 
have a good sensitivity if the hyperfine coupling lies in two 
regimes, see the rectangle I and the ellipse II. Inset (top left): 
contour plot for Ds versus a and B, in the case of k = 0.5/xs -1 
as above. Inset (bottom right): Ds, achieved by the optimal 
a, as a function of the radical pair lifetime r = k" 1 (in /is). 



assume that cf> — 0. The Hamiltonian then simplifies to 
H = B ■ Si + B ■ S-2 + aS^Iz- The nuclear spin polariza- 
tion remains unchanged during the spin dynamics. This 
is also true if the transverse component of the hyperfine 
coupling is small and thereby the spin flip is prohibited 
due to the energy mismatch. The nuclear spin can then 
be treated as inducing an effective magnetic field (de- 
pendent on its state) for the electron spin. If the nuclear 
spin is in the up (down) state |f) z (|l) 2 ), the effective 
magnetic filed is \az (— \az). Thus, we can written the 
effective Hamiltonian as H e = B ■ Si + B ■ S2 + bS^ 
where b = ±^a. After some calculations, we obtain 
the singlet yield as a function of 9 and a as $5(^,6) = 
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B) + § (1 + c) g(Bi - B), where c = cos(6» - 9'), g(x) = 
k 2 /(k 2 +x 2 ), B\ = (Bcos9 + b) 2 + B 2 sin 2 9, sm9' = 
Bshi9/B 1 , cos6»' = (Bcos9 + b)/B\. Therefore, the to- 
tal singlet yield for the magnetic with an angle 9 is 



<M0,f) + <M0,-f) 
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This analytic result for the singlet yield enables us to 
have a better understanding when a chemical compass 
can (not) have a good magnetic sensitivity. For instance, 
if the recombination rate k is too large, the radical pair 
will enter a quantum-Zeno regime and hence remain in 
the singlet state due to suppression of singlet-triplet tran- 
sition [15] . As a consequence, the singlet yield $s(9) — 1 
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is almost independent on the value of angle 9. A good 
sensitivity requires a moderate recombination rate to give 
enough time for magnetic field effects to occur without 
being severely degraded by spin decoherence. 

As verified by our numerical optimization results, we 
find that a chemical compass can work with a high and 
robust sensitivity if the hyperfine coupling lies in one 
of two regimes, see Fig.l. In the first regime (I), the 
hyperfine coupling is relatively large, i.e. a 3> B. In this 
case, we have the singlet yield as 

$ s (9) ~ i [1 + cos 2 {9 - 6')] ~ i (1 + cos 2 9) . (3) 

Therefore, the maximum and minimum singlet yields are 
$j max ' = $5(0,0) = 5, which corresponds to the coher- 
ent mixing between the singlet state \S) and the triplet 

state \T ) = ^(|t+} + |+t», and <^ min) = * s (f,a) = f 
corresponding to the mixing between the singlet state 
and all three triplet states. Thus, the magnetic sensitiv- 
ity is Ds = 2- This is exactly the sensitivity that can 
be achieved by a chemical compass which has one of the 
nuclear spin with a large dominant anisotropic coupling 
along a certain direction. For natural molecules which do 
not satisfy such conditions of hyperfine couplings, one can 
use magnetic nano-particles to design a metal-chemical 
hybrid compass and achieve the above sensitivity |25j . 
For very weak magnetic fields, the magnetic sensitivity 
is D S {B) = B 2 / [4 (k 2 + B 2 )] . This explicitly shows the 
dependence of the magnetic sensitivity on the the recom- 
bination rate k (i.e. the radical pair lifetime). Let us 
point out that the other benefit of having a large local 
effective field is to make the functioning of a chemical 
compass more robust against decoherence. In the sec- 
ond regime (see II in Fig.l) , the hyperfine coupling is 
smaller than the magnetic field, and the interplay be- 
tween them leads to an even better sensitivity (Ds = 0.4 
for k — 0.5/is" 1 and a — \B). The maximum singlet 
yield is obtained at 6 = j, where the hyperfine cou- 
pling is not large enough to cause the transitions from 
the singlet state to the triplet states (as a < B); while 
the minimum singlet yield happen at 9 = that corre- 
sponds to the coherent mixing between the singlet state 
\S) and the triplet state |T ). Let us stress that an opti- 
mally designed chemical compass can achieve a sensitiv- 
ity more than 10% even for a magnetic field that is one 
order magnitude weaker than the geomagnetic field (with 
the radical pair lifetime 2/is) - see inset in Fig.l. In fact, 
we find that a chemical compass may work at a magnetic 
field down to a few tens of nanotelsa if the radical pair 
life time is reasonably longer (e.g. ~ 10/is). This may be 
of practical interest if one wants to find applications for 
an artificial chemical compass. Finally, the robustness of 
such optimal configurations is shown in the contour-plot 
inset of Fig.l, i.e. the best performance is obtained in a 
relatively large range of the values of a and B. 
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FIG. 2. (Color online) (a) The singlet yield ® S (Q) as a func- 
tion of the angle 6 without (black square) and with (red circle) 
quantum control. Inset: simple control magnetic field C(t) (in 
mT) versus time (fis). A much weaker and shorter step func- 
tion (green area) for C(t) allows us to achieve a very similar 
sensitivity (~ 53%) as the optimal one (~ 54%), further show- 
ing the robustness of such external control, (b) The singlet 
ratio f s (t) for 8 = (green) and 9 = | (red), versus time (in 
us), with (dashed) and without (solid) quantum control. 

Extend magnetic sensitivity limit with simple quantum 
control. — As an application of quantum control in the 
context of radical pair reactions, it was shown that quan- 
tum dynamical coupling, which is a very specific kind 
of quantum control, can be used to disturb the func- 
tion of a chemical compass. The essence is to decouple 
the hyperfine coupling and thus malfunction the compass 
[T5] . The role of this special quantum control is destruc- 
tive. In contrast, here we show that by simple quan- 
tum control we can enhance the sensitivity of a chem- 
ical compass and go beyond the sensitivity limit with 
the optimal hyperfine coupling found above. We intro- 
duce an external control magnetic field (in the direction 
close to the x axis) that is a simple function of time as 
C(t) = J2k Ak sin (wjfct) +Bk cos (uikt) under some appro- 
priate amplitude and spectral constraints, which depend 
on the experimental conditions. Similar techniques have 
been recently proposed theoretically to prepare generic 
coherent or incoherent initial states of light harvesting 
systems |31j . Our numerical optimizations, based on var- 
ious algorithms as Nelder-Mead simplex direct search and 
subplex methods for unconstrained problems, shows that 
quantum control can improve the magnetic sensitivity 
significantly - see Fig. 2(a). For a chemical compass with 
the optimal hyperfine coupling, the sensitivity Ds (for 
the radical pair lifetime fc _1 = 2/is) is improved from 0.4 
to 0.54, which is even better than the sensitivity limit 
0.5 (that can be obtained only in the long radical pair 
lifetime limit without quantum control). For simplicity, 
we use two frequencies of oscillating fields with moderate 
amplitudes in order to generate realistic control fields 
as in the inset of Fig. 2 (a). The result is quite surpris- 
ing, as for other magnetic sensors the sensitivity tends to 
be degraded if one applies an additional magnetic field 
which may overwhelm the directional information of the 
magnetic field that we want to detect. However, in the 
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current case, the appropriate applied field (on both spins, 
and this makes it essentially different from the method 
with gradient fields) can instead amplify the directional 
effect of the weak magnetic field on the singlet yield. The 
phenomenon can be understood by looking at the evolu- 
tion of the singlet yield ratio over time for 8 = 0, see 
Fig. 2(b). In the first 2 ps, the control field is very small, 
the singlet state IS*) coherently transfers to the triplet 
state \Tq) , and then the increasing control field suppresses 
the back-conversion from \T ) to \S) (by inducing rapid 
transitions among the three triplet states) which reduces 
the minimum singlet yield. In the mean time, the con- 
trol field will keep the singlet state at 6 = leading to 
the maximum singlet yield. The above essential physics 
suggests that we can use other forms of control field (de- 
pendent on specific experimental conditions, which may 
not take the form of oscillating fields) to achieve a similar 
and robust sensitivity improvement. Indeed, motivated 
by this intuition, we have also applied a much shorter 
and weaker step-like magnetic field achieving essentially 
the same sensitivity - see inset Fig. 2(a). We remark that 
optimal control techniques can be exploited in a more 
general chemical compass with several nuclear spins to 
achieve a very high magnetic sensitivity enhancement, 
e.g. 60% for the FADH'-0 2 ~ molecule pi E3]. 

The optimization approach we use here can be also 
translated to design experiments with open-loop feedback 
control in order to achieve different targets, e.g. con- 
trolling the chemical yields and the radical pair lifetime. 
Once we gain more knowledge about different specific 
radical pair reactions from this type of spin chemistry 
experiments, it will be interesting to apply the obtained 
optimal control results to experiments on the actual bio- 
logical organisms, also verifying whether they really con- 
tain such investigated molecules. 

Decoherence- enhanced magnetic sensitivity of a chemi- 
cal compass. — The most prominent feature of avian mag- 
netoreception is that it can work at room temperature 
when various kinds of noise may exist. Under such con- 
ditions, quantum systems (here radical pair spins) suffer 
from decoherence and quantum features will be weaken. 
One will naturally expect that quantum decoherence will 
be unfavorable for the functioning of a chemical com- 
pass. Indeed, it was shown that if the decoherence rate 
is larger than the recombination rate, the magnetic sensi- 
tivity of a model compass will be dramatically degraded 
by general noises, while for a particular dephasing model 
(which is different from the model we study here), the 
singlet yield itself is entirely unaffected [2]. Here, we 
show that a chemical compass with an optimal hyper- 
fine coupling will not only have a good sensitivity, but 
is also quite robust against pure dephasing, and surpris- 
ingly, stronger dephasing may also enhance its magnetic 
sensitivity. For simplicity, we model noise as pure lo- 
cal dephasing described by the following Lindblad oper- 




FIG. 3. (Color online) (a) The magnetic sensitivity of the 
singlet, Ds, as a function of the dephasing rate 7 (ps -1 ). 
The radical pair recombination rate is again k = 0.5/is -1 . It 
shows the effect of uncorrelated (d = 0), partially correlated 
(d = 0.8), perfectly correlated (d = 1), and anti-correlated 
(d = — 1) dephasing on Ds for an optimally designed chemical 
compass, (b) Singlet yield &s(@) f° r uncorrelated dephasing 
with rates 7 = 0, 0.5, 2/is -1 , respectively. 
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and Lo = 



■ o~z ' , where 7 is the dephasing rate 
is Pauli operator. This quantum master equation 
approach is equivalent and can be directly mapped to a 
Liouville equation, see e.g. [33]. The parameter d char- 
acterizes how correlated the dephasing is, i.e. d = for 
uncorrelated dephasing and d = I for perfectly correlated 
one. We note that, since the distance between two radi- 
cals is of the order of several nanometers, the local energy 
fluctuations and thus dephasing may well be correlated. 
In Fig. 3, we plot the magnetic sensitivity as a function 
of the dephasing rate 7. It can be seen from Fig. 3(a) 
that the performance of such a chemical compass is quite 
robust and even enhanced by the presence of correlated 
dephasing. For the uncorrelated dephasing, the sensitiv- 
ity will first be degraded but will be enhanced again as 
the dephasing rate further increases. This result comes 
from the interplay between the coherent spin Hamilto- 
nian and the noise, which is evidenced in Fig. 3(b). A 
similar phenomenon happens for the anti-correlated de- 
phasing. Let us remark that, if the dephasing rate is 
extremely large and overwhelms all the coherent times 
scales, the spin dynamics is simply incoherent mixing of 
IS*) and |T ), and the sensitivity will eventually be de- 
graded as expected. 

Summary. — We investigate quantum limits for the 
magnetic sensitivity of a chemical compass by optimiz- 
ing its hyperfine coupling. We find that a model chemical 
compass with only one nuclear spin equipped by an opti- 
mal hyperfine coupling is already likely to give us a very 
high magnetic sensitivity. For such an optimally designed 
chemical compass, we demonstrate that its sensitivity can 
be further improved with simple quantum control. More- 
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over, its performance can be very robust and even better 
if the noise is correlated. In the case of uncorrelated 
dephasing, the magnetic sensitivity does not monotoni- 
cally decrease with the dephasing rate. These phenom- 
ena stem from the intricate interplay between quantum 
coherence and unavoidable noise. Therefore, our results 
may help us to understand what quantum coherent fea- 
tures of radical pair mechanism lead to chemical compass 
sensitivity, and provide useful design principles for a bi- 
ologically inspired artificial chemical compass. 
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